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An analytical method for getting new complex Hadamard matrices by using mutually unbiased 
bases and a nonlinear doubling formula is provided. The method is illustrated with the n = 4 case 
that leads to a rich family of eight-dimensional Hadamard matrices that depend on five arbitrary 
phases. 



O ■ 1. INTRODUCTION 

(N ; 

, Mutually unbiased bases (MUBs) constitute a basic concept of quantum information and play an essential role in 

■ quantum tomography, quantum cryptography, construction of dense coding schemes, teleportation, classical signal 
[ processing, etc. Its origin is in the Schwinger paper p3 [. 

^\ ' Two orthonormal bases in C^, A = (ai, . . . , a,]) and B = [hi, . . . , hd)^ are called MUBs if the moduli of the scalar 
product of theirs vectors ai and hj satisfy the relations |(ai,&j)| = 1 < i,j < d. An equivalent definition is that 

■ the product AB*o{ the two complex Hadamard matrices generated by A and B is again a Hadamard matrix, vifhere 
r*^ * denotes the Hermitian conjugate. 

It is well known that MUBs exist for any dimension d, a known result being that for any d > 2 there exist at least 
three MUBs, see [ij], [13l,0- 1^ '^^^ ^l^o established that the maximum cardinality of a set of MUBs in is less or 
equal to d + 1, with equality for d — , a power of a prime number p6l |. (isj ]. Q, [l7l |. For non prime numbers the 
results are very scarce such that even the simplest case, d = 6, is not yet solved, 

The technique for getting MUBs for p prime was given by Schwinger, [l^], who made use of the properties of the 
Heisenberg-Weyl group. Explicit constructions in dimensions one to five can be found in many papers, see e.g. [26j . 
0, i. 

^ • We are interested in MUBs that depend on arbitrary phases because we want to use them for the construction 
' of new complex Hadamard matrices in dimension d — To our knowledge the first who obtained such MUBs in 
CO dimension 4 was Zauner [l^l • Similar results are also given in [s] . 

The paper is organized as follows. In Section 2 we present some doubling formulas starting with those found by 
, Sylvester, [231] . and the nonlinear form of our formula, In Section 3 we show how to construct new eight-dimensional 
CNJ ' complex Hadamard matrices that depend on five independent parameters by using our nonlinear formula. 
C3 ' In Section 4 we discuss the Haagerup form of the equivalence of complex Hadamard matrices, [T(i|, and show that 
. behind it stand the orthogonality relations that lead to polygons in the complex plane. Because the unitary matrices 
[ ' are a particular case of normal operators we propose their characteristic polynomial as the true equivalence principle. 
I> . In Section 5 we collect our results written in the standard form introduced in We close the paper with a few 
concluding remarks. 



2. ORTHOGONALITY AND DOUBLING FORMULAS 



The orthogonality is an important concept that appears in many mathematically, physically or engineering contexts, 
and it is strongly related to the angle between two objects. The simplest case is that of a vector space where the 
angle is found from the usual inner product between two vectors, and if this product is zero the vectors are said to 
be orthogonal. This notion was extended to matrices, orthogonal designs, etc, and it is a strong tool used in many 
applications. One important problem is the use of an orthogonal object for the construction of a similar one whose 
dimension is bigger than that of the object one started with. It seems that one of the first who used this procedure 
was Sylvester in his paper, ,23], for doubling the size of Sylvester- Hadamard matrices. 

An important observation in p3| is the following: the known theorems relating to the form of the products of two 
sums of 2, or of 4, or of 8 squares must exhibit instances of orthogonal matrices of this nature. In this context he 
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observed that the entry wise product of the two-dimensional matrices 
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leads to the orthogonal matrix 
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The first matrix entering ([T]) is the simplest circulant matrix, and the second one is the core of almost all constructions 
involving Sylvester-Hadamard (real) matrices. In the four-dimensional case the three matrices have the form 
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(3) 



The first matrix entering ([3]) is also circulant being generated by the abelian group G = x 7^2, and the last matrix 
is the Williamson matrix, both of them being discovered avant la lettre. 

The parameters a, 6, c, d could be square (—1, 1) matrices of order m such that for every pair X,Y ^ (a, 5, c, d) 
the product XY*^ is symmetric, i.e. XY^ = YX* — (XF*)*, where t means transpose. If they satisfy the relation 
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cc 



dd = 4m/„ 



where /,„ is the m-dimensional unit matrix, the above W matrix is Hadamard of order 4to, [2^ 
starting point for construction of real Sylvester-Hadamard matrices by using Williamson matrix. 
For completitude we give also his solution for n = 8. 
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This is the usual 
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A little different formula for doubling the size of a unitary matrix was found by us in paper 
given three d x d unitary matrices. A, B, C, the 2d x 2d matrix 



(5) 



which says that 
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(6) 



is also unitary. Until now we used this formula only in the particular cases B — A and C = A for the construction of 
complex Hadamard matrices, i.e. we uncovered other two elementary block matrices similar to ([2]) 
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The doubling formula, i.e. the analog of Williamson matrix has a similar form 
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and so on. 

Another important point is that we can merge the above two constructions, ^ and (O, and obtain "hybrid" 
matrices of the form 
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which are orthogonal if A, B, C, D are circulant, or in a more general case if all XY* , X, Y G {A, B, C, D) are 
self-adjoint matrices, i.e. XY* — YX* = (XY*)*. The above three matrices could be also useful for finding real 
Sylvester-Hadamard matrices. 

In the following we do not explore these possibilities but we make use of formula ^ for obtaining new complex 
Hadamard matrices that could depend on more independent parameters than those involved in A, B, and C matrices. 
Here we restrict to a particular case, that of A, B, C coming from a set of MUBs, because it leads to a rich family of 
complex Hadamard matrices with a beautiful and unexpected symmetry. The first non trivial and interesting case is 
obtained when the three matrices have dimension four such that in this paper we discuss this case. However we start 
with the d = 2 case to see how the method works. 



3. NEW EIGHT-DIMENSIONAL COMPLEX HADAMARD MATRICES 



In the d ^ 2 case there is only one Sylvester-Hadamard matrix, up to equivalence, however we use this equivalence 
for taking a slightly different form for A, B, C matrices. Our choice was 
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where in the following i = a/— T. The above matrices have been multiplied by different two-dimensional diagonal 
phase matrices, and in this way each entry of CA*B matrix contains four terms. The resulting matrix gets unitary 
by choosing a certain relationship between the arbitrary phases, and this one can be done in four different ways, and 
all lead, up to a factor 1/2, to a matrix equivalent to the known result 
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By the parametrization change, a — !■ a + tt/2, the matrix ( 11 ip gets equivalent to the first complex Hadamard matrix 
that depends on an arbitrary phase, result obtained in |11| . 

In our case, d = 4, we use the three mutually unbiased matrices from (27j . in different combinations, and also the 
above result pT|) . These three matrices depend on arbitrary (real) phases and have the form 
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E^{w,z)=E^(x)-^E2[v,z) 
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where w = y — x \x\ the original definition, see |27j . but in the following we consider it an independent parameter. Our 
first choice is to take A = Ei, B = E2 and C = E3 and, after their multiplication at left, and/or right by diagonal 
phase matrices, each entry of CA*B matrix contains 2** = 16 terms. Because A, B, C matrices are Hadamard, the 
first problem to solve is the finding of those linear combinations of arbitrary phases such that the CA*B matrix should 
also be a complex Hadamard matrix. Each entry containing sixteen terms, the choice of independent phases is not 
unique, such that we can obtain many nonequivalent matrices. After that we have to check (impose) the Hadamard 
property for the resulting V matrix 

In this paper all the eight-dimensional matrices are defined up to the numerical factor 1/2^/2. If in the case = 2 
only one constraint between arbitrary phases was sufficient, in the d — A case the procedure implied at least four 
constraints, and the results are as follows: 
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D^^}{a,b,c,dJ) = 
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where the upper index, (5), means the number of independent phases, its lower part 8A means its dimension, and 
the first such a matrix. The above matrix has a beautiful symmetry that can be seen by looking at the last row and 
column and at the 2x2 inner sub-matrix 

" -1 -1 " 

-1 -1 

such that both these properties can be considered as a fingerprint of such matrices. We propose that the above form 
should be the standard form of a jacket matrix; more information about jacket matrices can be found in paper [l8!|. 

By using the symmetry of such matrices we can bring them to the above form, (|13p . with a plus sign in front of 
the first three exponential functions from the second column and row, and we consider it a manner to present this 
symmetry. However there are some cases when the used procedure leads to matrices whose second row has a little 
different form, and such a matrix is the following 
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(5) . . (5) 

Looking carefully at -Dg^ matrix one can observe that it can be brought to the same form as -Dg^ matrix if we permute 
between themselves the fourth and fifth columns, followed by a permutation of the second and the third rows, and 
respectively of the sixth and seventh ones. If on this new form we make simultaneously the changes a — >■ &, and & — > a, 
we find the matrix 
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Although the above matrix was obtained by using the usual definition of matrix equivalence for complex Hadamard 
matrices, jlOl] . -Dg^ and i'gp are not equivalent. The columns(rows) permutation, together with "benign" re- 
parametrization of the form, a b, and b a, lead to new non equivalent matrices. 

By using equivalent matrices to Ei, E2, E3 for A, B, C, obtained by permutations of rows/columns one get similar 
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matrices, that are not equivalent to ([T3 
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Another matrix with a httle different form is 
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Now we introduce into our game the matrix H^la), Eq. ()lip . that has a special dependence on the imaginary unit i 
by taking one or more of the A, B, C matrices, equivalent to i?4, each one depending on a different phase. Using the 
same procedure as above we got 
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It is easily seen that by the re-parametrization c — > c — 7r/2 and d ^ d — 7r/2, and, respectively, f ^ f — 7r/2, the 
matrices (dH) and ([20]) transforms into the matrix and respectively ([2^ . The matrices with, and without i, are 
not equivalent under the usual form of equivalence, however we will put them in the same class, as we will show latter. 
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Matrices with the same fingerprint as those found by us appeared in other different contexts. For example a similar 
matrix can be recovered from the complex Hadamard matrix, Pg, found in Q, that originally depends on three 
non-null complex parameters, and by permutation of rows and columns it can be brought to the form 
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that has the same symmetry. Thus we can recover from it a matrix which depends on five arbitrary phases, namely 
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and its transpose for f — a and d = coincides with the matrix (|2ip . 

Another example comes from a paper by Horadam, (l^ . who found a matrix, K4{i), that is equivalent to the back 
circulant matrix derived from a quadriphase perfect sequence of length 8, and such sequences are rare. 

By permuting rows and columns it can be written as 
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The above matrix has a similar symmetry to the previous matrices, and, by consequence, we can try to find the 
"natural" matrix it comes from. After few trials we found that the matrix (|23p is a particular case of the following 
complex Hadamard matrix that depends on five phases 
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If one makes the phases change b ^ b — tt/2, c c — 7r/2, d ^ d — 7r/2, f f — tt/2 one gets 
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The above matrices are new, see e.g. |24| and [J], where the only known result parametrized by five phases is the 
complex Hadamard matrix stemming from the Fourier matrix F^. 

All these matrices generate real and complex Hadamard matrices that depend on ±1 and ±i when the parameters 
(a, b, c, d, /) takes 0, ±7r/2 or tt values. There are 32 such matrices and the number of i entering each one is quantized 

q = 8n, n = 0, 1, 2, 3, 4. If we consider the matrix D^J(a, b, c, d, /), then there is one real matrix for n = 0, obtained 
when all the five parameters take 7r/2 value, eight matrices when n = 1, fourteen for n = 2, eight for n — i, and again 
one when n = 4. 

In the same time all the matrices generate 31 complex Hadamard matrices that depend on five arbitrary phases by 
doing re-parametrization of the form a — > a + 7r/2, where a denotes any phase entering a given matrix. 

One easily obtain similar matrices depending on Jour independent phases, the simplest method being to make zero 
a linear combination with integer coefficients of the five independent phases. However we used the same method by 
choosing other combinations of arbitrary phases entering the matrix C* AB and a few results are the following 
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The most interesting cases seem to be 
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because they cannot be brought by re-parametrization to matrices similar to (j26p and (j27p . This shows the importance 
of the Hi{a) matrix which allows to obtain more nonequivalent matrices that depend on four parameters. Thus the 
nonlinear term CA*B, whose each entry depends on 16 terms (phases), leads to non intersecting orbits depending on 
five, four or three phases, such that there do exist Hadamard matrices that depend on four or three phases that cannot 
be obtained from the generic five-dependent ones, e.g. fg^, by collapsing some parameters, or by re-parametrization 
of the form a — > a± 7r/2, where a denotes any parameter entering a given matrix. The above four matrices dep ending 
on four phases are new. In the known catalogue, there is only one such matrix, that found in paper jl9l |. which 
has another symmetry, showing that the set of all eight-dimensional complex Hadamard matrices has a rich structure, 
and the results obtained till now, including these ones, represent only partial results. 
An other example with three parameters is the following 
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4. EQUIVALENCE OF COMPLEX HADAMARD MATRICES 

The present day used form for the equivalence of complex Hadamard matrices is essentially the same as that for real 
matrices, and can be stated as follows. Two Hadamard matrices Hi and H2 are called equivalent, written Hi w H2, 
if there exist two diagonal unitary matrices Di and D2 and permutation matrices Pi and P2 such that 



Hi = D1P1H2P2D2 



(31) 



The above relation allows to bring any Hadamard matrix two its Sylvester form, '23|, i.e. to a matrix whose all entries 
from the first row and column are equal to unity. 

For complex matrices Haagerup suggested, [l3|, as a practical method for detecting (non)equivalence of two 
Hadamard matrices in more complicated cases, the use of the set of numbers 



{HijHkiHiiH, 



kj I «,J 



kA 



1, 



(32) 



for both the matrices Hi and H2^ where bar denote complex conjugation and Hij the entries of the matrix iJ, 
because the above product is invariant under equivalence pip . If the two sets of phases p2p coincide the matrices are 
equivalent. In some sense p2p is a local property that does not take into account the global properties of a complex 
Hadamard matrix. Let see what is behind it. 

For the first time such a product appeared in the paper [T^] , see also [l^ , which deals with the structure of the 3x3 
quark mass matrices in the standard electroweak model. These matrices for up and down quarks are diagonalized by 
unitary matrices U and U' and their product V = U U'* defines the so called Kobayashi-Maskawa unitary matrix. 
Jarlskog showed that the generic products Ira{VijVkiViiVkj) are invariant phases of the quark mixing matrix. In fact 
the invariant phases are directly related to the angles of the so called unitarity triangles, see e.g. @. 

In the general case the unitarity polygons are generated by the rows and columns orthogonality. For example the 
orthogonality relation of j and / columns is written as 



Ha = 



(33) 
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and if we consider its entries as vectors in a two-dimensional space the above relation, (j33p . defines a polygon in the 
complex plane. Because for an arbitrary n x n unitary matrix all the terms entering p3p could have different moduli 
the physicists divided the above relation by one term, let say HkjHki, such that one leg length is unity, and (|33|) has 
the new form 



1+ E 



HijHil 



= 



kj-nkl 



(34) 



For Hadamard matrices all legs length is unity if we redefine them without the factor 1/y/n, and then the generic 
term entering the above sum is HijHkiHkjHu which coincides with (j32p . Thus for complex Hadamard matrices the 
orthogonal relations (|33|) generate many polygons in the complex plane whose maximal dimension is n. All these 
polygons coming from all orthogonality relations define the fingerprint of a given complex Hadamard matrix, and we 
could say that two complex Hadamard matrices are equivalent if their fingerprints coincide. This fingerprint gives a 
global description of matrix properties. To see how this equivalence recipe works let us consider the matrix -ff4(a) for 
a = 0. 



1 1 f 1 

1 i -i -I 
1 -i i -I 
f -f -f 1 



(35) 



ft is easily seen that the fingerprint of the above matrix consists in the point 0(0, 0) of the coordinate axes, and the 
square OABC defined by ^ = (1,0), B = (1,1), C = (0,1). When a 7^ the Hi{a) fingerprint is the same point 
O, and a rhomb whose precise shape depends on the phase a. In conclusion we have a continuum of independent 
Hadamard matrices, all of them having the classical Hadamard form [ll| because the factor i can be absorbed in the 
phase a. 

We consider now two six-dimensional complex Hadamard matrices 
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1-1-1 1 i -i 
1 -1 -i -1 1 i 
1 1 -1 -i -1 i 
1 i 1 -1 -1 -i 
1 —i i i —i — 1 
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(36) 



According to relation p2p the above two matrices must be equivalent since both of them have the same number of 
i, namely 12, and the same number of —1, namely 9. However they cannot be equivalent because Di is a symmetric 
matrix, while D2 is self-adjoint. The relations p3p lead to similar fingerprints for both the matrices which consists 
in the preceding square OABC and the rectangle OABiCi with Bi = (1, 2), Ci = (0, 2), and the rectangle OA1B2C 
defined by Ai = (0, 2), B2 = (2, 1), and it is difficult to see if Z3g^ and Dqb are (non)equivalent matrices. In fact we 
have to compare all the 6 x 5 = 30 polygons generated by Dqa orthogonality, and, respectively, by Dqb- 

In the case of matrices D^J (a,b,c,d, f), * = A,...,L, the polygons can start with the point O, ending with 

octagons. For example the matrices D^g(a,b,c,d, f) and D^^,{a,b,c^d^ f) are not equivalent because the octagons 
defined by the orthogonality of the first and third columns have different shapes, etc. However, like the case of i?4(a) 
matrix, all matrices whose some entries contain products of the imaginary unit with arbitrary phases, and in case the 
imaginary unit can be absorbed by re-parametrization, will be considered as defining one continuous class. 

Defining the fingerprint class by all the orthogonality relations of the form ([55]) has an advantage over the present 
day used form because it implies that the matrices H, and H are all equivalent by reducing the number of non 
equivalent matrices. 

In the mathematical literature there is an other type of matrix equivalence for the class of normal matrices. A 
matrix N is normal if it commutes with its adjoint N* ^ i.e., NN* = N*N, for which the equivalence has a simple form. 
Every normal matrix is similar to a diagonal matrix D, which means that there exists a unitary matrix U such that 
N = U DU* , [111, p. 357, or [l^. Two important classes of normal matrices are the unitary and self-adjoint matrices 
such that we can say that two matrices, Ai and A2, are equivalent if and only if they have the same spectrum, or 
equivalently, the characteristic polynomials satisfy the relation, det{XIn — Ai) = det{\In — ^2), up to a multiplicative 
numerical factor, where det means the determinant of the corresponding matrix. 

For example the spectra of -DeA and Dqb matrices are as follows 



SpiD^A) = (^-l,l,iz(-\/6±iy30),i(-y6±zy3oj 
SpiD^B) = (-1,-1,-1,1,1,1) 



(37) 
(38) 
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such that the above matrices are not equivalent. 

The circulant matrix generated by the row {l,id, ~d, — z, —l/d,i/d), where d is one solution of the equation d^ — 
(1 — -\/3)rf +1 = 0, see Q, or the matrix (2.26) in [lOl, has the same spectrum as Dqb, so they are equivalent. 

In general the spectrum of a matrix cannot be explicitly found, however its characteristic polynomial can be, so 
the (non)equivalence can be established in most the cases. 

The characteristic polynomials of our matrices, up to numerical factors, are as follows 

Pn^^) (A) = (1 - A2)2 (2 A* - n/2 e'^A^ + ^2 e'^-^+^+'^+'^+^'^A - 2 e'(''+=+'i+/)^ (39) 
If one makes the re-parametrization d — > d + tt, which is equivalent to the interchange of the fourth and five columns 

(5) 

of the -Dg^ matrix, one finds 

P„(5!)(A) = (1-A2)2 (2 A^-v^e^'^A^- \/2e*(-''+''+=+'^+/)A + 2 (40) 



thus one can say that the corresponding matrices D^J and -Dg^^ are not equivalent. The same form for the relation 
(|40p can be obtained by using a re-parametrization with any other parameter, 6, c, or /, e.g., 6 — > 6 + tt, and all the 
four such obtained matrices are equivalent, as it easily seen from relations (p9)) and pO)) . 
For the next four matrices the polynomials are 

P^(5)(A) ^ (l-A2)2(2A'*-V2e*(^-"+^)A3-V2e^('^+'=+'')A + 2e*(°+^+''+/)) (41) 

P^(,^(A) = (1-A2)2 (2A4- V2e*'^A3 + V2e*(''+'^+-^)A-2e'(°+''+^'+-^)) (42) 

P^(5) (A) = (1 - XY (2A^ - %/2e"A3 + y2e*(-"+''+"+'*+/)A - 2 6*^''+=+''+^)) (43) 

P^(5)(A) = (l-A2)2f2A4-V2e"A3 + V2e'(^+''+-^'A~2e'("+'=+'^+-^'') (44) 

(5) 

The following matrix, D^p, leads to a more complicated polynomial 

P„(5)(A) = (1 - A2)2 (4 A'' - V2(2e*" + e'(^-'^+/) - e:f)X^ + 2(e'(''-'*) - l)(e*('=+'^) + e*(''+^))A2 

-V2 e''=(e*("+'') - e^("+'') - 2 e^(^+^))A - 4e*("+''+"+-'^') (45) 

The next matrix, P'gg, is interesting because by the re-parametrization c — > c — tt/2, and d — > c— 7r/2 it is transformed 
(5) 

into the matrix Dg^ . Its polynomial has the form 

P„(5, (A) = (1 - A2)2 (2 A" - v^e'^X^ - 72 e'(-'^+^+'=+'^+-^U + 2 e'(''+=+^'+-'')) (46) 



The formulas ([59]) and differ only by the sign of the last two terms, and could show us which are the matrices 
that are equivalent. In fact there are six such transformations involving any pair from the set (6, c, d, /) and all these 
four parameters, that preserve the form (|46p . If we make the transformation a ^ a + n in matrix (jl9p one gets 

P^,5i, (A) = (1 - A2)2 (2 A-* + v^e'^X^ + V2 e'(-"+''+'=+'^+/)A + 2 e*(''+=+''+^)) (47) 

The next two matrices lead to 

P^(5)(A) = (l-A2)2(e*(^+'i)-A2)(2A2 + \/2(ze*(-"+''+^')-e*")A-2ie*(''+-^)) (48) 

P^(5)(A) = (l-A2)2(e*(=+'')-A2)(2A2 + V2(e^(-'^+''+^'-e'")A-2e'(^+-^)) (49) 
whose spectra can be computed explicitly. A similar case is 

P^(5) (A) = (1 - A2)2(e^(=+/) + A^) (2 A^ - y2(e'(-°+''+'') + e")A + 2 (50) 

P^(,)(A) = (l-A2)2(e*(=+/)-A2)(2A2 + V2(e*(-''+''+'^)-e"')A-2e'(^+'*)) (51) 
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(5) 

Matrices of the form D^p seem to be rare, another example being the foUowing 
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and its polynomial is 



(1 - A2)2 (4 A'' + V2 [e'^b+'^-f^ - 2 e*"^ + e''')A^ - 2(e*("+'') + e*('=+-^'))(e*('*~-^) + 1)A2 



(53) 



which is similar to that of the o'fp matrix, see (pS)) . 

Looking carefully at the above polynomials we conclude that the Dg^, -Dg^ and D^q matrices generate a single 
class whose representative polynomial can be taken (j47|) since the polynomials ([39l) , (|40| , (|43l) , (|46p and (|47l) depend 

f d + /. The I?g^ and fi^'' matrices, 



5 + c + + / and 6 



and respectively 1507 and of}., define each one a new class, while Dan, o'f),, Dap, and £'i^r'* define five classes. 



on the same set of three distinct parameters, a, 

.(5) 

/gj ciiiix ' ^"^^i "i^^ ii^" vviiii^ ^8B' -^SC ^8E' ^8F """^ -^8L ^ 

Thus our proposal for equivalence of complex Hadamard matrices is the characteristic polynomial, which includes 
all the possible forms it could take by using simple re-parametrization, and/or row and column permutations. The 
standard form of these eight classes is given in the next Section. 

The spectra of matrices that depend on four phases are given by the following polynomials 

(54) 
(55) 
(56) 
(57) 
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which are similar to the preceding ones and show that the four matrices are not equivalent. 



("3") 

The matrix _Dg^ leads to 



P„(3,(A) = (1-A^)^ A-- 



1 
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(58) 



5. STANDARD FORM 



For writing the complex Hadamard matrices in their standard form, [2^, we need the real Hadamard matrices 

(5) . . (5) 

obtained by taking all the parameters zero in D^J , and the corresponding phase matrix -Rg, . There are four real 

Hadamard matrices 
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The corresponding RgJ matrices have the form 
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where 9 means zero. The standard form is 

ia,b,c,dJ)^h,o EXP (i ■Ri%,b,c,dJ)y Dfl (a, &, c, d, /) 



(66) 



(67) 



(68) 



= /i2 o EXP 



{i ■ iT!^^(a,6,c, dJ) 



(69) 



4'J(a, b, c, d, /) = /i3 o EXP z • R^i^ia, b, c, dJ)), D^^j>{a, b, c, d, /) = /^i o EXP » • i?^'J(a, b, c, d, /) (70) 



(a, 6, c, d, /) = /ii o EXP • (a, b, c, d, /)) , 4/ ^ c, d, /) = /14 o EXP • 4/ (a, ^ c, d, /) 



(71) 



i?fi(a, 6, c, d, /) = /i4 o EXP ( z • i?^'^^(a, 6, c, d, /) ) , i^^'^^ (a, 6, c, d, /) = o EXP ( * • (a, b, c, d, f) ) (72) 



?(5) , 



(5), 



In relation (I69p o denotes the entry wise Hadamard product. It is well known that matrices under Hadamard 
product and usual addition generate a commutative algebra such that all the usual functions like the exponential (p^ 
are well defined. 



CONCLUSION 



The use of the nonlinear doubling formula ^ led to an entire new class of eight dimensional complex Hadamard 
matrices which depend on five arbitrary phases. We did not make an exhaustive search for finding all such matrices, 
however we found many matrices, and tried to solve the problem of their equivalence by using the intrinsic properties 
of unitary matrices. By using the usual equivalence, pip , their number grows very quickly. Only the eight matrices 
(I69l) - (|72l) . by using re-parametrization of the form a — > a±7r/2 for the free parameters entering a given matrix, provide 
8 X 32 = 246 non equivalent matrices, that is a big number. 
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The use of orthogonality relations can solve a part of the equivalence problem, and an important result is that 
one matrix, its transpose and complex conjugate are equivalent. In this case comparing the fingerprints of two 
matrices is not an easy problem because the number of generated polygons if of the order n{n — 1) for each of the 
two matrices of order n. On the other hand the characteristic polynomial requires only one computation and, as the 
relations ((39| - (|53)) show, the characteristic polynomial depends only on a few combinations of the same set of five 
parameters (a, b, c, d, /), which are different for each class as defined before. This new equivalence definition includes 
the orthogonality relations equivalence since it is easily seen that because of the properties of normal matrices, the 
matrices, H, and H*, have the same characteristic polynomial, and that of H* is obtained by changing A — > 1/A in 
the previous ones, which leads to the complex conjugate spectrum. In this way the number of non equivalent matrices 
is significantly lowered. 

In principle the same equivalence could be applied to real Hadamard matrices. However a simple spectral compu- 
tation of the hi, /i2, hs, /14 matrices shows that only the matrices hi and h^ are equivalent, and hi is not equivalent 
to /i2 and /14, nor /12 is equivalent to /14. Their spectra are 

. (-,^,,^,l(l±£l^,^yi(L^^.^^'j ,73, 

where power means multiplicity. But taking into account that for real Hadamard matrices with dimension d = 
2, 4, 8, 12 there is only one matrix under the usual equivalence, we do not suggest the use of the new equivalence 
because it will cause dramatic changes in the field. 

In contradistinction to the real case the parametrization of complex Hadamard matrices is still at its infancy stage, 
such that the use of the proposed equivalence does not imply major changes. 

One important problem which deserves to be solved is the finding of all the possible classes of eight-dimensional 
complex Hadamard matrices that depend on five arbitrary phases. We used an ad-hoc method to find the above 
results, and presumably there are more. Closed related ones would be the description of the corresponding orbits, 
are they related, or not, to the form of characteristic polynomials, and solving the mystery why the use of the Hi{a) 
matrix leads to matrices which depend on four arbitrary phases, similar to p8| and (|29p . that cannot be brought to 
forms which do not explicitly depend on the imaginary unit. 
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